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1. Motivation 

The QCD phase diagram is conjectured to have a rich phase structure. At low temperatures, 
QCD has a vacuum and nuclear matter phase; at high temperatures and/or densities, QCD matter 
develops a qualitatively different phase where quarks are liberated from confinement - the so- 
called quark gluon plasma (QGP). While there is strong evidence for a crossover transition from 
the hadronic phase to the QGP for zero baryon chemical potential /Tg, there is no evidence for a true 
phase transition at higher densities. Lattice studies of QCD have aimed to extend the simulations 
to finite quark chemical potential /r = but the available methods are limited to /i/T < 1 due 
to the sign problem: Monte Carlo simulations sample a probability distribution and hence rely on 
the condition that the statistical weights are positive. In the conventional approach to lattice QCD 
based on the fermion determinant, the weight for the fermion determinant becomes complex as 
soon as the chemical potential is non-zero. The sign problem (more precisely in this context: the 
complex phase problem) is severe, prohibiting direct simulations for /r > 0 - which is also due to 
the fact that Monte Carlo is performed on the colored gauge fields. 

However, fhere is a represenfafion of laffice QCD which does nol suffer severely from fhe 
sign problem: in fhis represenfafion, fhe laffice degrees of freedom are color singlefs. The com¬ 
plex phase problem is reduced fo a mild sign problem induced by geomefry-dependenf signs of 
fermionic world lines. Such a “dual” represenfafion of laffice QCD has been derived for sfaggered 
fermions in fhe sfrong coupling limif, fhaf is in fhe limif of infinite gauge coupling g —)• o® [7]. In fhis 
limif, only fhe fermionic action confribufes fo fhe pafh infegral, whereas fhe action describing gluon 
propagafion is neglected. QCD af sfrong coupling has been sfudied exfensively since 30 years, bofh 
wifh mean field mefhods [1,2, 3,4, 5, 6] and by Monte Carlo simulations [7, 8,9, 10,11, 12]. Those 
sfudies have been limifed fo fhe sfrong coupling limif, which corresponds fo rafher coarse laffices. 
However, recenfly [13] we were able fo include fhe leading order gauge correcfions fo fhe partition 
function. The effecfs of fhese gauge correcfions on fhe phase diagram will be discussed below. 


2. The chiral and nuclear transition in the strong coupling limit 

The pafh infegral of sfaggered fermions in fhe sfrong coupling limif can be rewriffen exacfly as 
a parfifion funcfion of a monomer-i-dimer-i-flux system. The reformulafion proceeds in fwo steps: 
firsf fhe gauge links (gluons) are infegrafed ouf, which confines fhe quark fields vr(x) info color sin¬ 
glefs, fhe hadrons: fhose are fhe mesons M{x) = t//(x)t/r(x) and fhe baryons 

B{x) = ^£ip2i^\|/i^{x)\|/h\|/i^(x). In fhe second step, also fhe quarks are integrated ouf, which al¬ 
lows fo express fhe parfifion funcfion via infeger variables: 

Zsc{mq,^l)= ^ (2-1) 

b={x,p) ,c n^- ^ i 

V-—V--( 

meson hoppings condensate M, baryon hoppings 


wifhkfo G {0,... 3}, G {0,... 3}, G {0, ±1}. Since fhe quark fields are freafed as anli-commuling 
Grassmann variables in fhe pafh infegral, fhe infegrafion realizes a Pauli exclusion principle called 
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Figure 1: The Phase diagram in the strong coupling limit (left), as measured in a Monte Carlo 
simulation, compared to the standard expectation of the continuum QCD phase diagram fright). 
Both diagrams are for massless quarks. 

the Grassmann constraint: 

nx+ ^ (ky{x) + ^\£y{x)\\ =3. (2.2) 

v=±o,...,±i ■' 

This constraint restricts the number of admissible configurations {kh,nx,£} in Eq. (2.1) such that 
mesonic degrees of freedom always add up to 3 and baryons form self-avoiding loops not in contact 
with the mesons. The weight w(£,/f) and sign o{i) = ±1 for an oriented baryonic loop I depend 
on the loop geometry. The partition function Eq. (2.1) describes effectively only one quark flavor, 
which however corresponds to four flavors in the continuum (see Sec. 4). It is valid for any quark 
mass. We will however restrict here to the theoretically most interesting case of massless quarks, 
lUq = 0. In fact, in this representation the chiral limit is very cheap to study via Monte Carlo, 
in contrast to conventional determinant-based lattice QCD where the chiral limit is prohibitively 
expensive. 

Eor staggered fermions in the strong coupling limit, there is a remnant of the chiral symmetry 
1 / 55 ( 1 ) C SUiiNf) X SUR{Nf). This symmetry is spontaneously broken at T = 0 and is restored 
at some critical temperature Tc with the chiral condensate (yw) being the order parameter of this 
transition. As shown in Eig. 1 (left), we find that this transition is of second order. This is analogous 
to the standard expectation in continuum QCD with Nf = 2 massless quarks, where the transition is 
also believed to be of second order. Moreover, both for our numeric finding at strong coupling and 
for the expectation in the continuum, the transition turns into first order as the chemical potential is 
increased. Thus the first order line ends in a tricritical point, which is the massless analogue of the 
chiral critical endpoint sought for in heavy ion collisions. 

In fact, at strong coupling, the zero temperature nuclear transition at /ig c — >^3 is intimately 
connected to the chiral transition, and they coincide as long as the transition is first order. The 
reason for this is the saturation on the lattice due to the Pauli principle: in the nuclear matter 
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Figure 2: Example of an ^(j8) diagram. On an excited plaquette, color singlets can also be com¬ 
posed of quark-quark-gluon or antiquark-gluon combinations. Whereas in the strong coupling limit 
baryons are pointlike, they become extended objects due to the gauge corrections. 


phase at r = 0, the lattice is completely filled with baryons, leaving no space for a non-zero chiral 
condensate to form (in terms of the dual variables, there is no space for monomers on the lattice). 
This is certainly a lattice artifact which disappears in the continuum limit, where the nuclear phase 
behaves like a liquid rather than a crystal. 

The ultimate question is whether the tricritical point at strong coupling is related to the hypo¬ 
thetical tricritical point in continuum massless QCD. If we can establish such a connection numer¬ 
ically, this would be strong evidence for the existence of a chiral critical endpoint in the pt-T phase 
diagram of QCD. To answer this question, it is necessary to go away from the strong coupling limit 
and incorporate the gauge corrections, which will lower the lattice spacing and eventually allow to 
make contact to the continuum. 

3. Gauge Corrections to the strong coupling phase diagram 

Lattice QCD in the strong coupling limit is dehned by the the lattice coupling j3 = ^ 0 as 

g —)• oo. Going away from the strong coupling limit is realized by making use of strong coupling 
expansions in j3. We have recently shown how to incorporate the leading order gauge corrections 
^(j3) [13]. In a nutshell, the strategy is to compute link integrals at the boundary of “excited” 
plaquettes, which correspond to gluonic excitations. Introducing a vaiiable qp G {0,1} to mark 
the ’’excited” plaquettes P, the & (jS) partition function can be expressed in a similar fashion as 
Eq. (2.1) with modihed weights vv (for details see [13]): 

{nA,i,qp} X b i P VO/ 

We can sample this partition function by the same algorithm (variant of the worm algorithm) as for 
=0, adding a Metropolis accept/reject step to update the plaquette variables qp. These simulations 
have been caiTied out for Nt = 4 and various lattice volumes Na =4,6,8,12,16 to perform hnite 
size scaling and to measure the phase boundary as a function of the chemical potential. In contrast 
to the strong coupling limit, where the color singlets are entirely composed of quarks and anti¬ 
quarks, including the gauge corrections allows color singlets to be composed of quark-quark-gluon 
or antiquark-gluon color singlet states, as shown in Eig. 2. There are two qualitatively new features 
that arise when incorporating the corrections: 
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Figure 3: Phase boundary in the fi-T plane extended to finite j8. The backplane corresponds to the 
strong coupling limit j3 = 0. The second order phase boundary is lowered by increasing j3. We 
do not observe a shift of the chiral tricritical point. However, the nuclear critical endpoint (CEP), 
determined from the baryon density, moves down along the first order line (extrapolated to T = 0 
to guide the eye) as j3 is increased. 

1. Baryons are point-like in the strong coupling limit, the lattice spacing is too coarse to re¬ 
solve the internal structure of the baryon. Including the gauge correction, baryons become 
extended objects, spread over one lattice spacing. 

2. The nuclear potential in the strong coupling limit is of entropic nature, where two static 
baryons interact merely by the modification of the pion bath. With the leading order gauge 
correction, pion exchange is possible as the Grassmann constraint is relaxed: on excited 
plaquettes, the degrees of freedom in Eq. (2.2) add up to 4 instead of 3. 

These features will have an impact on the phase boundary. In Eig. 3, the effect of the gauge 
corrections is shown. We find fhaf fhe second order phase boundary is lowered, as expecfed because 
fhe critical femperafure in laffice unifs drops as fhe laffice spacing is decreased wifh increasing j3. 
However, we find fhe chiral fricrifical poinf and fhe firsl order fransifion fo be invarianf under fhe 
^(j8) correcfions. We wanf fo sfress fhaf fhere are acfually fwo end poinfs, which splif due fo 
fhe gauge corrections: fhe second order end poinf of fhe nuclear liquid-gas fransifion is fraced by 
looking af fhe nuclear densify as an order parameter. We expecf fhe nuclear and fhe chiral firsf 
order fransifion fo splif, such fhaf af T = 0 fhere are fhree differenf phases insfead of fwo phases (as 
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shown in Fig. 1 right). The nuclear phase is in the continuum distinct from the chirally restored 
phase. As a first evidence for this splitting, we find that the nuclear critical end point separates 
from the chiral tricritical point. 

4. Relation between the strong coupling phase diagram and continuum QCD 

In Fig. 4 we speculate how the separation of the first order transitions could he realized at 
larger values of jS. Moreover, we can distinguish at least three scenarios (A,B,C) on how the chiral 
tricritical point depends on jS. These scenarios start from the same phase diagram in the strong 
coupling limit, but have different continuum limits at j3 —)■ oo (a —)■ 0). In all three scenarios, a 
tricritical point exists at /i = 0, j3 > 0: it must exist because the finite-temperature /i = 0 transition, 
which is of second order for j3 = 0, is of first order for j3 = oo, following the argument of [14] 
which applies to the continuum, Nf = 4 theory. 

1. In scenario (A) the chiral transition remains first order for all values of /ig. Hence the tricrit¬ 
ical line turns towards /i = 0 at some finite ■ 

2. In scenario (B) the chiral transition weakens and hence turns second order, but strengthens 
again to turns first order at larger jUg. 

3. In scenario (C) the chiral transition weakens and remains second order. In that case the 
tricritical line bends towards larger /r and eventually vanishes at some finite 

In order to discuss the relation between the phase diagram in the fJt-T plane for Nf = 4 massless 
quarks with the more physical scenario Nf = 2 + 1 with 2 massless up and down quarks and one 
physical strange quark, we show phase diagrams in the Nf-jJ. plane. Interpolating between integer 
numbers of massless flavors Nf and Nf + 1 can be realized by decreasing the mass of an additional 
flavor from infinity to zero. In all scenarios it is assumed that for Nf = 2, the chiral transition is 
second order, and that there is a tricritical strange quark mass separating it from the Nf = 3 
first order transition, as shown in the so-called Columbia plot. Fig. 5. Note that whether Nf = 2 is 
indeed second order and thus whether exists and also whether it is larger or smaller than the 
physical strange quark mass is still under debate [15]. The standard scenario of QCD in the chiral 
limit, as shown in Fig. 1 (right), corresponds to scenario (B) in Fig 4. However, the non-standard 
scenario (C) is supported by Monte Carlo simulations for imaginary chemical potential and analytic 
continuation [15, 16]: these studies suggest (at least for small chemical potential) that the chiral 
transition weakens with chemical potential, making the Nf = 3 first order region in Fig. 5 to shrink 
with increasing jUg. This should also be the case for Nf = 4. 

A last comment on staggered fermions is in order: one of the lattice artifacts is due to the 
way this discretization solves the so-called fermion doubling problem: At strong coupling, there 
is effectively only one quark flavor, whereas in the continuum limit the same action describes 
4 flavors due to the fermion doubling. Instead of 15 Goldstone bosons that are present in the 
Nf = 4 continuum theory, there is only one Goldstone boson at strong coupling, since the others 14 
receive masses from lattice artifacts (called taste-splitting). In the determinant-based approach, the 
problem is solved by “rooting”: taking the root of the fermion determinant to reduce the number of 
flavors from 4 to 2 (and the number of Goldstone bosons from 15 to 3). This strategy is not available 
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(A) 


(B) 


(C) 



Figure 4: Top row: Various scenarios for extending the phase diagram in the strong coupling limit 
(j8 = 0) toward the the continuum limit (j8 —?■ oo). All three scenarios assume that the nuclear and 
chiral transition split, and that at /i = 0 the chiral transition is of first order (since in the continuum 
Nf = 4). In the strong coupling limit, the chiral transition at /r = 0 is second order (corresponding 
to Nf = 1 and the doublers decoupled), hence there must be a tricritical point at some It 

is an open question whether the tricritical point at strong coupling is connected to this tricritical 
point at (Jeft), or connected to the speculated tricritical point in the continuum (center) or 

terminates at some finite (5 at T = 0 (right). 

Bottom row: the corresponding scenarios for the finite temperature chiral transition in the p — Nf 
phase diagram, showing the possible relation of the tricritical point at Nf = 4 with those at Nf = 
2 + 1, assuming the chiral limit for the light quarks and a physical strange quark mass. The p-Nf 
is limited by the line Pc{T = 0), beyond which chiral symmetry is restored. Left: For Nf = 4, the 
transition is of first order for all values of p. Center: The tricritical point at V/ = 4 is is connected 
to the tricritical point at V/ = 2 + 1. This would be evidence for the existence of the critical end 
point in the QCD phase diagram for physical quark masses. Right: The Nf = 4 first order region 
does not extend to Nf = 2 + 1, where it remains second order. This second order transition turns 
into a crossover immediately as mu,ni^ > 0, so in this scenario there is no chiral critical end point 
at physical quark masses. 

Vf=0 

Figure 5: The Columbia plot with the assumption 
which implies that the chiral transition is second order for 
Nf = 2. The arrow points towards the V/ = 2 + 1 chiral light 
quark masses and physical strange quark mass as denoted in the 
bottom row of Fig. 4 in between Nf = 2 and Nf = 3. 
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in our dual-variable approach. Although the strong coupling limit has effectively only one flavor, 
the residual chiral symmetry is that of a Ay = 4 continuum theory, with one true Goldstone boson 
which even persists when the chiral anomaly Ua (1) is present for j3 > 0. This is in contrast to a 
genuine Nf = I theory in the continuum which has no Goldstone bosons at all. The chiral anomaly 
breaks the chiral symmetry explicitly, driving the chiral transition into a crossover (corresponding 
to the lower right comer of the Columbia plot Fig. 5). Hence the deconfinement transition at 
Nf = 0 is most likely completely separate from the chiral transition for Nf > 2, as shown in all 
three scenarios Fig. 4 (bottom). 


5. Outlook for future investigations 


There are various ways to discretize fermions on the lattice, with staggered fermions and 
Wilson fermions the most widely used for thermodynamics studies. They describe the same physics 
in the continuum limit only. At finite lattice spacing, and in particular at strong coupling, both 
discretizations are quite different. In particular, the spin and the kinetic term of the fermion action 
are treated very differently. A dimer-i-flux representation is also possible for Wilson fermions. 
Such a representation was so far only determined for lattice QED [17, 18], since the Grassmann 
integration is much more involved for > 1. 

As a matter of principle, for both lattice discretizations, the gauge action can be incorporated 
order by order in j3. There are however technical difficulties that remain to be solved. A new 
strategy to study both lattice discretizations on a par is to expand both in systematically in j3 and 
the inverse quark mass by making use of a Hamiltonian formulation [19]. The partition function is 
then expressed by a Hamiltonian composed by operators: 


Z = Tr[/-^], 




2 T. 

^ fj) Qi 


4 = (4)’ (5 1) 


where the generalized quantum numbers Qi (spin, parity,flavor) are globally conserved, and near¬ 
est neighbor interactions are characterized by the operators ■^Qfx)'^Qiy)’ which raise the quantum 
number Qi at site x and lowers it at a neighboring site y (see [ 19] for the case of Nf = 1,2 for stag¬ 
gered fermions). For both staggered fermions and Wilson fermions, the matrices Jq contain vertex 
weights which are the crucial input to sample the corresponding partition function. The plan for 
the future is to do so with a quantum Monte Carlo algorithm. Comparing both fermion discretiza¬ 
tions order by order in the strong coupling expansion will help to discriminate lattice discretization 
errors from the genuine physics, in particular with respect to QCD at finite density. 
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